Spatially Anisotropic Heisenberg Kagome Antiferromagnet 
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We study the quasi-one-dimensional limit of the spin- 1/2 quantum Heisenberg antiferromagnet on the kagome 
lattice. The lattice is divided into antiferromagnetic spin-chains (exchange J) that are weakly coupled via inter- 
mediate "dangling" spins (exchange J'). Using one-dimensional bosonization, renormalization group methods, 
and current algebra techniques the ground state is determined in the limit J' ^ J. We find that the dangling 
spins and chain spins form a spiral with 0(1) and 0{J' /J) static moments, respectively, atop of which the chain 
spins exhibit a smaller 0[{J' / J)'^] antiferromagnetically ordered component along the axis perpendicular to the 
spiral plane. 



I. INTRODUCTION 

The nearest-neighbor Heisenberg antiferromagnet on the 
kagome lattice, a two-dimensional network of corner-sharing 
triangles, is one of the most geometrically frustrated mag- 
nets. Frustration suppresses the magnetic ordering tendency 
and leads to an extensive classical ground state degeneracy. 
Order-by-disorder effects lift the degeneracies in the classical 
system and are believed to select the coplanar \/3 x \/3 pattern 
as the ground stateJ-^iM 

The spin- 1/2 quantum kagome antiferromagnet is much 
less understood. Exact diagonalizatio n^i^i^'^'^- '°' ' ' and series 
expansion studie s'^''-^ indicate the absence of long-range mag- 
netic order, but the precise nature of the ground state remains 
mysterious. For small systems, numerical simulations reveal 
a large number of singlet states below a small (and possi- 
bly vanishing) spin gapi ^'^-'"^ This observation has lead to the 
speculation that the ground state of the kagome lattice might 
be a gapless critical spin liqui d'^''^''^ or a particular type of 
valence-bond crystal that exhibits many low-energy singlet 
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(see also Refs.l21landl22 



state! 

Recently, two new candidate materials for an ideal spin-1/2 
kagome antiferromagnet have attracted considerable atten- 
tion. First, the mineral ZnCu3(OH)6Cl2, also known as Her- 
bertsmithite, realizes structurally undistorted, magnetically 
isolated kagome layers with Cu^+ moments on the lattice 

.jt^. 23.24.25.26.27.28.:- 



2^ Neither magnetic ordering nor spin freez- 



ing has been observed for this material down to the lowest cur- 
rently achievable temperature of 50 mK, which is well below 
the energy scale of the antiferromagnetic interaction. ^"^-^^'^^ 
While Herbertsmithite is structurally perfect, Dzyaloshinskii- 
Moriya (DM) interactions and a small number of impuri- 
ties might complicate the experimental study of the ideal 
quantum kagome system.^- Secondly, there is Volborthite 
Cu3V207(OH)2-2H20, a spin-1/2 (Cu2+) antiferromagnet, 
whose magnetic sublattice consists of well-separated kagome- 
like planes 4ii^2i^^i^ This material has a monoclinic distortion, 
which deforms the equilateral kagome triangles into isosceles 
triangles leading to a difference between two of the nearest- 
neighbor exchange constants (J') and the third one (J). Sim- 
ilarly as in Herbertsmithite, the spins do not order down to 
1.8 K, an energy scale fifty times smaller than the exchange 
coupling strength in Volborthite.^' However, at very low tem- 



peratures evidence for a spin freezing transition has been 
reportedi^ The anisotropy ratio of the exchange couplings, 
a = J/ J', could not be determined experimentally so far, 
but the differing side lengths of the kagome triangles seem to 
favor a > 1. A recent comparison of exact diagonalization 
calculations with thermodynamic measurements suggests that 
the spatial anisotropy of the exchange couplings is small, and 
that additional interactions beyond the nearest-neighbor cou- 
plings might be present in the mineral Volborthite.^^ 

In this paper, motivated by the renewed interest in the 
kagome systems and the recent experiments on Volborthite, 
we investigate the spatially anisotropic version of the quan- 
tum kagome antiferromagnet. We shall focus on the quasi- 
one-dimensional limit, J' <C J, where the model consists 
of quantum-critical spin-1/2 chains weakly coupled together 
via intermediate "dangling" spins (see Fig. [Til. In this situ- 
ation the competition between quantum fluctuations and the 
strong geometric frustration of the kagome lattice is particu- 
larly keen. The anisotropic quantum kagome antiferromag- 
net has been studied previously by a variety of techniques, 
all of which employ perturbation theories in some "artificial" 
small parameter Examples include large- A^ expansions of the 
Sp(A) symmetric generalization of the modelr^i^ a block- 
spin perturbation approach to the trimerized kagome lattice^ 
and semiclassical calculations in the limit of large spini^i^ 
Our approach is complementary to these studies in that it of- 
fers a fairly controlled analysis of the quasi-one-dimensional 
limit using powerful field-theoretical methods^2il2ili that have 
originally been developed for the investigation of quantum 
critical systems in one dimension. Recently, the quasi- 
one-dimensional version of the kagome antiferromagnet in a 
strong magnetic field, assuming that the intermediate spins are 
fully polarized, has been studied using similar techniques'^ 
The present paper treats the case of zero external field. 

Our approach rests on the assumption that in the quasi-one- 
dimensional limit, J' <C J, the intermediate spins order at a 
temperature scale Tg much higher than the ordering temper- 
ature Teh of the weakly coupled chains. This is justified a 
posteriori by our finding that the effective interaction among 
the interstitial spins, which sets the temperature scale Tg, is of 
order {,J'Y / J, whereas the most relevant effective interaction 
among the weakly coupled spin-1/2 chains, which determines 
Teh, is of order {J')*/J^. Consequently, we can divide the 
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FIG. 1: (Color online) Spatially anisotropic kagome lattice with 
nearest-neighbor exchange J (blue) among chain spins (S) and J' 
(black) among chain spins and interstitial spins (s). 



theoretical analysis into three separate stages. First, we de- 
rive the effective interaction among the intermediate dangling 
spins using perturbative renormalization group (RG) transfor- 
mations in the time direction. These considerations are com- 
plemented by numerical estimates of the induced short dis- 
tance couplings among the dangling spins. Second, we an- 
alyze the ground state of the resulting, spatially anisotropic 
triangular lattice of dangling spins as a function of first and 
further-neighbor interactions. We find that in the ground state 
the interstitial spins form a rotating spiral with a small (and 
possibly vanishing) wave vector parallel to the chain direction. 
Third we determine the most relevant interchain interactions 
using a symmetry analysis and RG considerations. Besides 
the effective interchain interaction, the spiral magnetic field 
produced by the intermediate spins induces another perturba- 
tion to the system of decoupled Heisenberg chains. Finally, 
we analyze these perturbations to the fixed point of the inde- 
pendent (decoupled) spin- 1/2 chains with the help of operator 
product expansions (OPEs). 

The ultimate result of our analysis is that all spins order, 
with the non-coplanar configuration shown in Fig. [21 in which 
the interstitial and chain spins predominantly form coplanar 
spirals with a wave vector {q, 0), but with a reduced 0{J' / J) 
static moment on the chains. The chain spins are weakly 
canted out of the plane, with the 0[( J'/ J)^] normal compo- 
nents forming an antiferromagnetically ordered pattern. The 
precise value of q cannot be reliably determined from our 
analysis, but we expect g <C 1 and q = is a distinct 
possibility. In the q = Q case, the state becomes coplanar. 
This ordered state differs from those found in two other re- 
cent studies^il^ using other methods, but there are similari- 
ties. These are discussed in Sec. lVIII 

At first glance, it might be counter-intuitive that the intersti- 
tial spins and the chain spins order in a nearly ferromagnetic 
fashion among themselves rather than in an antifeiTomagnetic 
pattern. But it has to be kept in mind that the ordering of the 
quasi-one-dimensional version of the kagome antiferromagnet 
is driven by the weakly coupled interstitial spins, which order 
at a larger energy scale than the chain spins. There is no a pri- 
ori reason that the effective interaction among the interstitial 
spins should be purely antiferromagnetic. In fact, it turns out 
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FIG. 2: (Color online) (a) Three-dimensional perspective view of 
the spin ordering pattern. The interstitial spins (red arrows) form a 
coplanar spiral with an 0(1) local moment in the x-y plane. The 
components of the chain spins (gray arrows) in the x-y plane also 
form a spiral but with static moment of 0\J' j J\ and antiparallel to 
the interstitial spins. The out-of-plane components (2 direction) of 
the chain spins are non-zero but even smaller - 0[(J'/J^)^] - and 
ordered antiferromagnetically along and between the chains, (b) Top 
view of the spin ordering pattern. The vertical component of the 
chain spins is indicated using + (upward pointing) and — (downward 
pointing). 



that there is an effective ferromagnetic interaction along the 
diagonal bonds connecting nearest neighbor interstitial spins 
(see Sec. llVb . Together with an effective antifeiTomagnetic 
interaction along the horizontal bonds connecting neighbor- 
ing interstitial spins this ultimately leads to the spiral order 
of the dangling spins. The spiral order of the chain spins can 
then be understood as arising from the linear response to the 
local field of the ordered interstitial moments. 

The remainder of the paper is organized as follows. Sec.HIl 
describes the lattice Hamiltonian, its symmetries, and its low- 
energy field theory description. In Sec. [Ill] we analyze the 
low-energy interactions using symmetry considerations and a 
perturbative renormalization treatment. In Sec.|IV]we use nu- 
merical methods to estimate the interaction strength among 
the interstitial spins and derive an effective model for the dan- 
gling spins on the triangular lattice. Sec. [V] deals with the 
ground state analysis of the triangular lattice. The perturba- 
tive analysis of weakly coupled chains in a spiral magnetic 
field with a small wave vector, 5 ^ 1, is presented in Sec. I VII 
and our conclusions and a discussion of the relation to other 
results are given in Sec. lVIII Some technical aspects of the RG 
calculations are relegated to Appendix [A] For completeness. 



we study in Appendix |B] the ordering of the weakly coupled 
chains in the presence of a spiral field with a large wave vec- 
tor. 



II. MODEL DEFINITION AND LOW-ENERGY 
HAMILTONIAN 

The Hamiltonian of the Heisenberg antiferromagnet on the 
anisotropic kagome lattice (see Fig. [Til is given hy H = 
Hq + V, where Hq describes the decoupled set of chains with 
nearest-neighbor antiferromagnetic Heisenberg interactions 



^0 — J 7 ^ Si^y ■ Si + l 



y 



(la) 



i,y 



and V is the interaction among the chains and the intermediate 
spins 



V = J'^ 



^2j±i,2yTi 



• {S2i,2y + S2i±l. 



2y 



+ S2i,2y^l + S2i±l,2y^lj ■ (lb) 



For the interstitial spins we shall use the symbol s, whereas 
the chain spins are denoted by S. The anisotropic kagome 
lattice has rotation, translation and reflection symmetries. The 
translational subgroup is generated by the translations Ti and 
T2, which move the lattice by two units along the horizontal, 
and one of the diagonal axes, respectively. The rotational sub- 
group consists of TT rotations about the lattice sites and about 
the centers of the hexagons. We can distinguish two types 
of reflection symmetries: reflections Ri about a vertical line 
through a midpoint of a chain bond (link parity), and reflec- 
tions i?2 about a horizontal line passing through the interme- 
diate dangling spins (see Fig.fTli. 

Let us now describe the chains in the continuum limit, 
which is applicable as long as J ^ J'. In this limit the low- 
energy properties of the antiferromagnetic spin- 1/2 chains 
are governed by the Wess-Zumino-Novikov-Witten (WZNW) 
SU(2)i theory, and the chain spin operator Si,y can be de- 
composed into its uniform My{x) and staggered Ny{x) spin 
magnetizations 



ao[My{x) + {-irNy{x)]. 



(2a) 



where x ~ iao and gq denotes the lattice spacing. The uni- 
form magnetization can be written in terms of left- and right- 
moving SU(2) currents. My = Jy,R + Jy.L- Another impor- 
tant operator describing low energy properties of the spin- 1/2 
chains is the staggered dimerization Sy [x), which is defined as 
the continuum limit of the scalar product of two neighboring 
spins 



S^+l^y -^ {~lTey{x). 



(2b) 



The scaling dimension of these continuum operators deter- 
mines the relevance of the operator in the RG sense with 
respect to the Luttinger liquid fixed point of the decoupled 



chains. The uniform magnetization My has scaling dimen- 
sion 1, whereas both the staggered spin magnetization Ny and 
the staggered dimerization Sy have scaling dimension 1/2. 
The microscopic lattice symmetries of H, Eq. ([T]i, leave an 
imprint on continuum description (|2]i. That is, the action of the 
space group symmetries on the continuum operators is given 
by 

Ti : M -* M, N ^ +N, e ^ +e, (3a) 

i?i : M ^ M, N ^ -N, e -^ +e, (3b) 

R20T2 : M ^M, AT -> -iV, s -^ -e, (3c) 

where, for brevity, we have suppressed the chain index. Other 
symmetry operations on the continuum fields are either trivial, 
or can be rewritten as a product of the above transformations. 
The three continuum fields M, N, and e form a closed op- 
erator algebra with respect to certain operator product expan- 
sions, which are widely used in the literature i^^-^"'"^''"^-^i^"^'"^^'"^^ 
For example, the right-moving SU(2) currents Jji satisfy the 
following chiral OPEs 



J^(a:,T)4(0) 



(5'^V(87r2) 



(uT — ix + aofTr)^ UT — ix + aoffr 

(4a) 



with imaginary time t, ar ~ sgn t, the short-distance cut-off 
ao, and spin velocity u = TTJao/2. Similar relations hold for 
the left-moving spin currents Jl . The product of J^ and TV 
can be expanded as 



JUx,t)N\0) 



47r(ur — ix + ooctt) 



(4b) 



The above equalities are understood to be valid only when in- 
serted into correlation functions and in the limit where the two 
points (.T, r) and (0, 0) are close together. These current alge- 
bra relations will allow us to compute one-loop RG equations 
by purely algebraic means (see Sees. Hill andlVIb. 

Using relation (l2al i. we can derive a naive continuum limit 
of the interaction among the chains and the intermediate spins, 
Eq. dTbli. First, we note that the intermediate spins Sjj^i 2!,_i 
couple symmetrically to S2i,2y and S2i+i,2y Hence, the stag- 
gered spin magnetization enters only via its first derivative in 
the continuum version of the interaction V, and we can set 

S2^.2y + S2^+l,2y ^ 2aoM2y(2x) - ^9, ATa^ (2x) . (5) 

In the above, and the following, we use the notation that the 
derivative in dxN2y{2x) is with respect to the full argument, 
i.e., more explicitly 

d,N2y{2x) = dxN2y{X)\^^^^ . (6) 

With this, the interaction Eq. dTbl i reads V^ = Vi + V2, with 

Vi = 71 51 *2.±i,2„T^ • iM2y{2x) + M2y^l{2x)] , 

(7) 



^,v 



V2 = 72 2^(±)s2x±i,2.yT^ • d- [^2,/(2x) + N2y^l{2x)] , 



x,y 



where the bare coupling constants are given by, 71 = 2 J'oq 
and 72 = —J'a^/2, which follows from substituting (|5]l 
into ( fTbl i. We have retained the next-to-leading interaction 
V2 as it will produce, in combination with the leading term 
Vi, relevant interchain couplings with respect to the fixed 
point of the decoupled chains (see Sec. HiH) . The contin- 
uum description of the interaction V , Eq. (|7]i, is necessarily 
invariant under the symmetry transformations, Eq. (O. In 
particular, we note that V does not contain a contribution 
S2x+i- 2y-i ' -^2i/(2x) which is forbidden by the symme- 
tries of the spatially anisotropic kagome lattice (Unk parity 
Ri sends N -^ -TV). 



III. RENORMALIZATION GROUP TREATMENT 

Here, we study the effective low-energy interactions among 
the dangling spins {H/\), as well as those between nearest 
neighbor chains (ych), using both a symmetry analysis and 
perturbative renormalization group transformations. The goal 
is to understand which interactions amongst interstitial spins, 
amongst chains, or betwixt the two, are most relevant, and 
further, at what energy scales they become important in de- 
termining the low-energy physics. Technically, the Hamilto- 
nian ([T]i with the chains treated in the continuum limit is for- 
mally rather similar to a Kondo lattice. As the Hamiltonian re- 
tains some local character due to the intermediate spins, in the 
framework of the RG approach, this problem develops only in 
time, and the energy is the only variable that is being rescaled 
as the RG progresses. Such an RG scheme is similar in spirit 
to the one employed in the context of a single ir npurity cou- 
pled to a Luttinger liquid (see, for example, Refs. [47B8B9h . 



A. Symmetry analysis 



dimerization Sy and the staggered magnetization Ny. There- 



fore, we find 



Vc 



ch 



y 



/ dX {jNNy ■ Ny + l + -f^EyEy + l} , 



(8) 



where the value of the coupling constants -fN and 7^ will have 
to be determined by microscopic calculations. Using Eqs. (O 
it is straightforward to check that these are the only terms with 
lowest possible scaling dimension 1 that satisfy the symmetry 
requirements of the spatially anisotropic kagome lattice. 

In addition to these most relevant interchain interactions, 
it is necessary to consider also a few less relevant terms, as 
these will arise at lower order in the renormalization group 
treatment below, and are important for generating the more 
relevant terms in Eq. (O above. These are 

Kh ' =J2 r^^ {^ONd^Ny ■ da:Ny+l + jAlMy ■ My+l} . 

y -^ 

(9) 
Next, we turn to interactions among the interchain spins 
i^A- Considering only terms that arise at 0[{J'Y], we find 
that H/s is given by Heisenberg interactions among dangling 
spins whose y coordinates differ by at most one unit. The cou- 
pling constants of these Heisenberg interactions are restricted 
by the symmetries of the lattice. With these conditions, i^A 
can be conveniently written in the form (see Fig.|3ll 



i^A - 2^ j2rS2x: 
a:,y,r>0 



±i,2./T| ' *2a:±i+2r,2yT3 



(10) 



+ 2^ ■^\2r+l\S2x±^,2yTk ' *2a;±|+2r,2y±i • 

Since the couplings J^ decrease in magnitude with increasing 
\r\, we can truncate the above sums over r after the first few 
terms. 



Before proceeding with the RG derivation of the low- 
energy effective interactions, we first write down the most 
general form of T4h and H^ that are allowed by the symme- 
tries of the spatially anisotropic kagome lattice. Such a gen- 
eral symmetry consideration will reveal all relevant symmetry 
allowed interactions that are expected to be generated through 
RG transformations, when all nominally irrelevant terms are 
taken into account. The space group symmetries needed for 
this analysis have been discussed in Sec.|II][see Eqs. (O]. 

We begin with the interchain Hamiltonian V^h- In principle, 
the number of allowed interchain interaction terms is infinite. 
However, only a handful of terms will be important. Most 
significant are those terms which are most relevant in the RG 
sense (with respect to the Luttinger liquid fixed point of the 
decoupled chains). This amounts to two-chain interactions 
involving no derivatives and only continuum fields that have 
small scaling dimension. Amongst these, we may further re- 
strict ourselves to nearest-neighbor chain interactions, as the 
magnitude of further-neighbor chain interactions is expected 
to decrease with separation distance. The continuum oper- 
ators with the smallest scaling dimension are the staggered 



B. Renormalization group results 

To determine the fluctuation-generated corrections to the 
low-energy effective action we perform a perturbative RG 
analysis of the interactions V , Eq. (|2]i, to one-loop order The 
perturbation theory is formulated by expanding the partition 



function Z 
plings 



/' 



-Sq-J drV 



up to quadratic order in the cou- 



Z c^ f e-'^" 1~ f drV + ^Tf dTidT2V{Ti)V{T2) 

(11) 

Here, 5*0 denotes the fixed point action, T is the time-ordering 
operator, and n is the imaginary time. Implicit in Eq. (fTTI) 
is a regularization, such that the integrals appearing in the 
expansion are taken only over the regions in which no two 
times are closer than some short time cut-off a = ao/u. The 
RG proceeds by increasing this cut-off infinitesimally from 



a to ba, where b 



cdt 



and di! > is the usual logarith- 



mic change of scale. To do this, all pairs of times ri , T2 such 



that a < |ri — T2I < ha must be fused using the operator 
product expansion, and the integral over ti — T2 in this range 
carried out. One thereby obtains a new partition function with 
renormalized interactions and the increased cut-off. We then 
perform an additional trivial rescaling of time and space coor- 
dinates for to restore the original cut-off: 



My{x,T) -> b^^My{x/b,T/b), 

Sx.yir) -^ S^,y{T/h). 



(12) 



We note that the spatial rescaling is not required to restore 
the cut-off but is natural and convenient for the continuum 
fields of the chains, which are described by a Lorentz invari- 
ant field theory. By contrast, we obviously cannot rescale the 
coordinates of the discrete spin operators of the interstitial 
spins. This difference leads to the appearance of an explicit 
RG length scale in the effective couplings between the chains 
and interstitial spins, Vi and V2, in the renormalized Hamilto- 
nian: 

Vl{t) = 71 ^ «2.±i,2,T| • [M2y{2x,) + M2y^^{2x,)] , 
x,y 

V2{i) = 72E(±)^2.±i,2,Tr^4^2,(2x,) (13) 

+N2y^ii2xi)\, 

where xg = a;e~^. 

With this, the derivation of the RG equations is straightfor- 
ward, and follows closely the methods used in Refs.|47|48|49l 
Details of the calculations can be found in Appendix|A] Tak- 
ing into account y = T^i + V2 to one loop, we obtain RG 
equations for both the interchain couplings and the couplings 
between the interstitial spins. The former are 
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Note that the strongly relevant interactions jpf, 7^ are not gen- 
erated at one loop from 71, 72. Instead, the former are gener- 
ated from the sub-dominant Jon,1m terms. This is the reason 
the latter terms needed to be included in our treatment. 

The RG equations for the interstitial spin couplings are 
functional, insofar as they describe the flow of the full set of 
interactions J7r. We find 



dj2r 

de 



J2r 



2 2 

-lMire-')^+lNire-')'^, (15a) 
u u 



where 



^"M - E 3^ 



f^^^ 47r2 [a + cr'i2r]^ 
lN{r) = 8Cn: 



a{a'^ -8r2) 



(q:2 + 4r2)5/2 



(15b) 
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FIG. 3: (Color online) Exchange paths of the triangular lattice 
formed by the intermediate spins s of the spatially anisotropic 
kagome lattice. Solid black lines denote the nearest-neighbor ferro- 
magnetic exchange ^1 (a - b and b - c bonds), dotted black lines the 
nearest-neighbor antiferromagnetic exchange >72 (a - c bond). The 
further neighbor interactions ^3 (a - d bond) and >74 (a - e bond) for 
the spin Sa are illustrated with red arrows. The chain bonds of the 
anisotropic kagome lattice are indicated in light blue (cf. Fig. [Hi. 



where Cat is the amplitude of the (N°-N°-) correlator Anal- 
ogous expressions can be derived for the interactions J7|2r+i|> 
i.e., the second term in H/s, Eq. ( fTOl i. 

In passing, we note that the form of the RG flows in 
Eqs. ( fT4l ) is highly constrained by the symmetry of the full 
(spatially anisotropic) kagome lattice, and this is necessary to 
avoid the generation of strongly relevant interactions at this or- 
der. By contrast, for the kagome strip of extension one in the 
y direction, there is no translational symmetry T2, and conse- 
quently symmetry ( l3cb is absent. This allows for the appear- 
ance of a term in the renormalized Hamiltonian proportional 
to / dx£y{x, r), which is indeed generated proportionally to 
7172 in that case (compare with Azaria et al., Ref. Is^, and 
see Ref. pli) . This term is strongly relevant and generates a 
gapped dimerized state in that model. 



C. Implications 



The RG flows in Eqs. (1141151) give considerable insight into 
the emergent energy scales in the problem. Though interac- 
tions between the chains and between the interstitial spins are 
both generated from the bare interactions 71 , 72, their charac- 
ters are distinctly different. Crucially, we see from Eqs. ( fT4l ) 
that amongst the chains, only marginal and irrelevant interac- 
tions - 7a jv, 7m - are generated at this order. Since the bare 
values 71 (£ = 0) '--^ 72(^ = 0) are both of 0[J'], this includes 
all effective generated interactions up to 0[(J')^]. A simple 
scaling argument shows that the low-energy excitations of the 
chains can be modified from those of decoupled chains only 
on energy scales parametrically smaller than {J'Y ■ 

To understand this, let us consider the solutions to Eqs. ( fT4l ) 
in somewhat more detail. The most important observation is 
that the strongly relevant interactions, 7^ , 7^ [third line of 
Eqs. (fT4bl. once they are generated, grow rapidly and expo- 



nentially with i. They are therefore dominated by the ef- 
fects of their "source" terms (proportional to 7ajv7A/) at small 
i, which are most amphfied under the growth. The source 
terms at small scale are themselves of order {J'Y / J, as de- 
scribed above. Therefore we expect that the relevant interac- 
tions will be, at least until they renormaUze to large values, of 
order ( J')^/ J^. Indeed, the ultimate effect of all this short- 
scale renormalization is completely equivalent to including 
fluctuation-generated "bare" interactions of this same order. 
This procedure has been described in considerable detail in 
prior publications, in which the formal manipulations closely 
parallel those used here. See Refs. 39 41 for further informa- 
tion. Therefore we will in Sec. [VT] simply take the relevant 
couplings to have the appropriate values: 



7w(£ - 1) = + 



'\4 



2(J0 
47r4 J3 ' 



7e(^ 



1 =- 



47r4j3' 



(16) 



By contrast, the interactions Jr generated in Eqs. dTsl l at 
0[( J')^] amongst the dangling spins will clearly modify their 
energetics at this same order, since they act within an oth- 
erwise completely degenerate manifold. More formally, this 
follows directly from the RG flows, Eq. ( fTsl l, which shows that 
the interstitial interactions will grow to 0(1) values (actually 
to any fixed value) by a scale Je^ ~ Jili^li) ^ i^'Y' / J ■ 
Therefore we see that, as claimed in the introduction, the or- 
dering (as we shall see, this is the result of these interactions) 
of the interstitial spins occurs at an energy scale at which the 
chains are unaffected by their couplings to the interstitial spins 
and each other 

Having understood the hierarchy of energy scales resulting 
from the RG, we can proceed to try to understand in more de- 
tail the ordering of the interstitial spins. In principle, one may 
do so by integrating the flow equations for the J,, couplings, 
Eqs. ( fTSl l, from £ = to the ordering scale, I = 21n(J/J'). 
This is equivalent to performing the time integrals not over 
the infinitesimal shell in the vicinity of the running cut-off, 
but instead from microscopic times a up to times of order 
a{J/ J'Y. Inspection of expressions in Eq. dlSbt shows that 
these integrals in fact are dominated by times of order the mi- 
croscopic time. Indeed the upper limit of the time integration 
can be extended to infinity without significant modification of 
the result. 

The dominance of short times has physical significance: 
it implies that the generated interstitial spin interactions are 
in fact controlled by high-energy physics of the Heisenberg 
chains, i.e., correlations of the chain spins on the scale of a 
few lattice spacings. These short distance correlations are not 
precisely captured by the continuum bosonization/current al- 
gebra approach. Therefore, if an accurate determination of 
these dangling spin interactions is desired, we must abandon 
(for the moment) the field theory methodology and return to a 
direct study of the lattice model. Conversely, because the gen- 
erated interactions are dominated by short distance physics, 
we expect that no serious divergences will arise in a micro- 
scopic calculation. In the following section we show how per- 
turbative and numerical methods may be used to perform the 
necessary calculations and obtain the fluctuation generated in- 
teractions amongst the interstitial spins. Using these precise 



results we determine the interstitial spin ordering in Sec. |V] 

Following this, in Sec.[VT]we will return to the lower energy 
physics resulting from the fluctuation-generated interactions 
7a jv, 7a/ in Eqs. (fl4] i. which will ultimately lead to ordering 
amongst the chain spins as well. 



IV. NUMERICAL ESTIMATES OF INTERACTIONS 
AMONG DANGLING SPINS 

The dynamical structure factor of the spin- 1/2 antiferro- 
magnetic Heisenberg chain, S{q, lo), can be related to the cou- 
pling constants of the effective interactions among the dan- 
gling spins, Jr, by means of second-order degenerate pertur- 
bation theory in the interaction V, Eq. ( fTbj i. Using available 
numerical data for S{q,uj), both in the two-spinon approxi- 
mation and for finite-size systems, this gives an estimate for 
the sign and relative strength of the interactions Jr. 

To formulate the perturbation theory in V we denote the 
ground state of the unperturbed Hamiltonian Hq, Eq. dTal i. by 
|0) and introduce the ground-state projection operator P = 
|0)(0|. With Q = 1 - P, any given state |*) of chain (S) 
and interchain (s) spins can be written as |$) = \'i'o) + I^q), 
where I^Pq) = P|*) and |*q) = Ql*). We note that I^Jfo) 
describes the state with the chain spins in the ground state, 
while the state of the s-spins remains arbitrary, i.e., |^o) = 
|0)|{s}). Following along the lines of Refs. i42ll52l we obtain 
an eigenvalue equation for |^o) 



Hn+PV 



1 



1-RQV 



i?F |*n)=S|*n 



(17) 



with the resolvent R ~ {E — Hq)^^. Eq. ( [TtI i is highly non- 
linear in the energy E, since the left-hand side of this equation 
depends on E through the resolvent R. Performing a perturba- 
tion expansion on the operator (1 — RQV)~^ gives at second 
order in J' 



Ho + PVRoV]\'<5Jo)=E\'<5/o 



(18) 

where Rq = (Eq — Hq)^^ is a function of the ground state 
energy Eq of the decoupled chains. Multiplying Eq. ( fTSl ) 
from the left by (0| and inserting the perturbation V, Eq. ( [Tbb . 
yields, after some algebra, a Schrodinger equation for the in- 
terchain spins alone 

H^\{s}) = iE-Eo)\{s}), (19a) 

with the effective interaction among interstitial spins 



x,y,r>0 

+2(J')^^^(|2r + l| 



'2a:±i,2j/^i ' ^2x±^+2r,2y^^ 



x^y^r 



*2a;±i,2y^i ' ^2x±^+2r,2y±i 

(19b) 
Here, A{r) is the linear combination 

Air) = GM(r - 1) + 2GM{r) + GA/(r + 1) (20) 



of the ground state expectation values 

GM{r) = (01^2°.,, F \ ^2.+r.,|0)- (21) 

In deriving equation ( fT9] ) we made use of the fact that the ex- 
pectation value ( |2TI ) is independent of both the x and y coordi- 
nates, due to translational invariance. The linear combination 
in Eq. ( l20l i arises, because every interstitial spin s is coupled 
to pairs of chain spins S. From inspection of Eq. (fT9l l we 
find the following expressions for the exchange couplings Jr, 
which have been defined in Eq. (fTOl i. 



J2r = A{J'fA{2r), 
J2r+i = 2{J'fA(2r + l), 



(22) 



with r > 0. The factor of 2 difference between the first and the 
second lines in Eq. ( l22l i accounts for the fact that interstitial 
spins s with the same y coordinate, say, i/ — i, are connected 
via both the y and y — 1 chains. 

Finally, by inserting a resolution of identity in Eq. ( |2TI) . 






1 



-Eq — E„ 



{n\Si,+r.M, 



with JTi) denoting the eigenstates of i/o, we realize that 
Gm{t), Eq. (ISTT i. is just a spectral representation of the zero- 
frequency Matsubara spin Green's function 



GM{r) = - / dT(YrS''{2x + r, r)5'^(2x, 0)). (23) 

The zero-frequency Green's function G]\i{r) in turn is 
connected to the dynamical structure factor S{q,uj) via a 
Kramers-Kronig transform 



du;' I \qS{q,oj')^^^. (24) 



GM{r) 



Inserting the above equation into the definition of A{r), 
Eg. (I20I 1. gives 



Air) 



^ Tdu;' rdqcos^1s{q,co'f-^.i25) 



Having related the interstitial exchange interactions J7r 
to the dynamical structure factor, we evaluate numerically 
S{q, w) in the subsequent sections to obtain precise estimates 
for the first few interstitial couplings Ji • • • J4. 



with the fundamental integral 

e* cosh(2t) cos(4pt) - 1 



Hp) 



dt- 



t 



coshisinh(2t) 



(27) 



■ sin q and 



and with the lower continuum boundary uj^ 

the upper continuum boundary lo]^ = tt sin |. The auxiliary 

variable p is a function of q and lo 



p{q,Lu) = — acosh - 

TT 



I UJ^{q) - Luliq) 



[(q) 



(28) 



Inserting 5*2(9, tj) into formula ( |25] |. we obtain A{r) in the 
two-spinon approximation 



A2{r) 



A r^ 

-^ dq cos^ - cos(qr) (29) 

T^ Jo 2 






Expression (I29b can be evaluated numerically in an efficient 
way, provided one splits off the singular part from inte- 
gral ( |27] | (see Ref. l53b . In doing so, we obtain for the ratio 
between the couplings J2 and Ji the following result within 
the two-spinon approximation 



J2 
Ji 



2^2(2) 
^2(1) 



-0.7446, 



(30a) 



and find that J'l is ferromagnetic, while J'2 is antiferromag- 
netic. Similarly, the magnitudes of further neighbor interac- 
tions are estimated to be 



J3 
Ji 



-0.1909, 



J4 
Ji 



-0.2312. 



(30b) 



with J3 < and J4 > 0. 

The two-spinon intensity accounts for about 73% of the to- 
tal structure factor intensity.^^^ The remaining part is carried 
by states with a higher number of spinous, and it is believed 
that the four-spinon intensity together with the two-spinon in- 
tensity cover about 98% of the spectral weight.^"* In principle, 
it would be possible to evaluate the ratios in Eq. ( l30b within 
the four-spinon approximation. However the involved numer- 
ical integrals are rather expensive to compute. Instead, we 
shall use finite size results for the structure factor to obtain a 
second estimate for the magnitude of the couplings Jr. 



A. Two-spinon dynamical structure factor 

First we compute the couplings Jr using the two-spinon 
approximation for S{q,u!). The two-spinon contribution to 
the dynamical structure factor can be explicitly written as (see 
Refs.pSlandK 



Miu^-uj], (26) 



1 g-I{p{l:'^)) 

S2{q,iLi) = , „ =6[( 



27r 



\MAm) 



B. Finite size results 

We compute A{r) using numerical data from the ABACUS 
database^^ for the dynamical structure factor in a finite system 
of A^ = 500 sites. These numerical data were obtained using a 
method based on the Bethe ansatz framework, which involves 
a summation over the so-called determinant representations 
for form factors of spin operators on the latticCi^- In these 
computations the momentum delta functions are smoothed out 



8 



by including a finite broadening 77 ^ 1/N. Using the numer- 
ical data with a smearing 77 = 0.01 we obtain for the ratio 
J2/J1 



Ji 



2Ajv(2) 
An{1) 



-0.7013, 



(31a) 



with J7i < and J'2 > 0, in accordance with Eq. (I30ab . The 
estimated magnitudes of further neighbor interactions are 



J3 
Ji 



^ ~ +0.2349, ^ 



Ja 
Ji 



-0.1453, 



(31b) 



where jTs is ferromagnetic and J74 is antiferromagnetic. As 
anticipated from the RG treatment in Sec. HUB I we observe 
that the magnitude of the coupling strengths decreases with 
spin separation distance. 

Comparing Eqs. dSOl l to Eqs. (ISTT i. we infer that the val- 
ues of J2I J\ agree well (within 10%), whereas the agree- 
ment between further neighbor interactions is worse. Since 
the two-spinon approximation misses 27% of the structure 
factor intensity, we suspect that the true values of the cou- 
pling strengths are closer to the finite size results (ISTT i than to 
the two-spinon estimates. In the next section, we are going to 
work out the interstitial spin ordering using results ( |30] l and 
(ISTT i and neglecting any longer ranged interactions with Jr>A- 



V. SPATIALLY ANISOTROPIC TRIANGULAR LATTICE 

As the discussions in the previous sections have shown, 
the ordering of the dangling spins s, which occurs at an en- 
ergy scale 0[(J')^], is independent of the chain ordering. At 
the energy scale 0\{J'Y\ the spins s, which form a spatially 
anisotropic triangular lattice^2iS2iS2iffi, are described by the ef- 
fective Hamiltonian i^A, Eq. ( fTOl i, with antiferromagnetic 
interactions along the horizontal bonds (J2) and ferromag- 
netic interactions along the diagonal bonds ( Ji), see Fig. |3] 
Neglecting third-nearest neighbor and longer ranged interac- 
tions, we can truncate the sums over r in Eq. ( fTOl i after the 
first two terms and obtain 

[y] («J> 



for 



(u)' 



where [ij] and \if\ denote the diagonal bonds and (ij) and 
(if) the horizontal bonds connecting first and second-nearest 
neighbors, respectively (see Fig. O. The values of the cou- 
pling strength J\- ■ ■ Ji are given by Eqs. (13013 It . 

This is a non-trivial spin-1/2 quantum model. However, 
it is also one which has been heavily studied, at least in the 
nearest-neighbor limit. In this case, over the vast majority of 
the phase diagram, the quantum ground state agrees with the 
classical one. This is particularly true when the couplings are 
such that the classical ground state is "close" to ferromagnetic 
(indeed the fully polarized ferromagnetic state is of course an 



exact eigenstate, as usual). We therefore expect that a classical 
analysis is reliable, and pursue it below. 

The classical phase diagram of Hamiltonian ( [32] i is found 
by replacing the spin operators with classical coplanar spiral 
vectors, 



Si = xcos{q ■ Ti) -|-ysin(q- n), 



(33) 



and minimizing the energy, which amounts to minimizing the 
Fourier transform of the exchange couplin g^''^^ 



J{q) 



hi 



Jij cos [q- (r, -Tj)], 



(34) 



where r^ denotes the position of site i in real space, and x and 
y are two orthogonal unit vectors. In our case, Eq. (|34l i gives 

J{q) = 2 Ji cos Qa; cos Qy + J2 cos{2qj;) (35) 

+2J3 cos{3qx) cos Qy + J4 cos(4(3'j;). 

with Ji, J3 < and J2, J4 > 0. 

For simplicity, let us first analyze the minima of Eq. 
for the case J'i~Ji = 0. This gives q = {q^ , 0) with 



Qx 



0, 



J2<\Jl\/2, 



arccos[-Ji/(2J2)], J2 < \Ji\/2. 



(36) 



That is, for J2 < |>7i|/2 the classical ground state is fer- 
romagnetic, whereas for J'2 > \Ji\/2 the ground state is a 
spiral state rotating along the x direction with the wave vec- 
tor q — {qx,0). [To study how quantum fluctuations alter 
the classical ground state in the case J3 — J4 ^ Q we have 
performed a linear spin-wave analysis (not shown). The mag- 
netization decreases smoothly from 1/2 to zero as J2/\Ji\ is 
increased, with a kink at J2/\Ji\ = 1/2.] 

In the case of non-zero further neighbor interactions, J^ ^ 
and J4 ^ 0, there is no simple explicit expression de- 
scribing the global minima of J{q), Eq. ( [35] ). Numerically 
we find that the minimum of J{q) with the coupling values 
given by Eq. dJTI ) (the finite size results) occurs at q = 0, 
i.e., the ferromagnetic state is the ground state. However, 
the coupling parameters of Eq. ( [3T1 ) are rather close to the 
boundary of a spiral phase in the coupling parameter space 
{Jr I f = 1, • ■ ■ 74}. In particular, the coupling parame- 
ters in Eq. ( [30l ) (the two-spinon result) yield a spiral state with 
q ~ 27r(0.08, 0). Therefore, we shall consider in what follows 
a cycloidal spiral ground state with wave vector q ~ {q^, 0), 
with qx small, which includes as a limiting case the ferromag- 
netic state [qx = 0), that is. 



So ( ices 



(37) 



O I Qx 

2qxX ± — 



+ ysin 



n , Qx 

2qxX ± — 



Here sq < 1/2 is the local static moment (staggered magneti- 
zation) of the spiral state. 



VI. COUPLED HEISENBERG CHAINS IN A SPIRAL 
MAGNETIC FIELD 

We are now in a position to address the spin-ordering pat- 
tern of the chains. We assume that the interstitial spins have 
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ordered into the spiral state given by Eq. ( l37b . and focus on the 
lower temperature scale of order ( J')^/ J'^, at which the chain 
spins become affected by relevant interchain interactions. The 
unperturbed system of antiferromagnetic Heisenberg chains is 
given by (note that we let a; ^ a; — 1/2 compared to SeclVTl 



Hq - J^S^_^^y •S'^+i,y. 



(38) 



The perturbations are: (i) coupling between chain and inter- 
chain spins, described by ( |40l i. (ii) marginal "backscattering" 
term, which is already present for a single Heisenberg chain, 
and (iii) the fluctuation-generated interchain interactions ([8]). 
Since the interchain spins form a two-dimensional ordered 
spiral state, described in the preceding section, the main effect 
of their interaction with the chain spins S is described by the 
spiral magnetic field vector {sx) given by (note that we let 
qx -^ q compared to SeclVT) 



SQ[xcos{qx) +yaui{qx)], 



(39) 



with q small. The spiral magnetic field introduces the follow- 
ing perturbations to this system of decoupled spins 



Hs = 



S = ^'X!{'^2a.-l/2,a-(s2:r-l+S2:.) (40) 

+ S2x + l/2,y ■ {S2x + S2x + l)\ 



x^V 



7 , ''•a:+l/2 ■ ^x + l 



with 



h,, 



/2,y, (41) 



2/io cos((7/2) [x cos{qx) + i) sm{qx)] , (42) 



and the field strength ho = sqJ'. 

The simplest effect of the spiral field, as for any field, is 
to induce a corresponding spiral magnetization in linear re- 
sponse. This implies that, generically, there will be static com- 
ponents of the chain spins in the spiral (x-y) plane. Since the 
spiral field is proportional to J', and the chain susceptibility 
is generally proportional to 1/ J, these static components are 
of order J' /J. 

This simple linear response, however, is not the complete 
story. A much more detailed analysis is needed to resolve the 
more subtle effects of the spiral field beyond linear response, 
in conjunction with the inter-chain couplings described in 
Sec. Uni The nature of this more complete analysis depends 
crucially upon the magnitude of q. When q is small, the spins 
respond in a way which is similar to the response to a uni- 
form magnetic field. A systematic approach is then possible, 
in which the Hamiltonian is transformed into the slowly ro- 
tating frame in which the external field is uniform. This is 
tractable because for small q the non-Heisenberg interactions 
induced by the change of frame are weak. Because of the ex- 
pected smallness of q, we focus on this case in the remainder 
of this section. 

In the opposite limit of a "large" wave vector, which is 
incommensurate with the dominant fluctuations of the ID 



Heisenberg chain, i.e. q, Jtt — g| ^ 0(1), the field weakly 
couples to the spin chain. Indeed, the leading 0{h) effects av- 
erage out over space, and instead only sub-dominant terms are 
induced at 0{h^). In this limit these weaker 0{h^) terms are 
crucial in determining the final state of the system. Though 
this limit is actually conceptually simpler than the opposite 
one, it is technically challenging, because the effects of spiral 
field are determined by very short-wavelength properties of 
the Heisenberg spins. As a result, we are not unambiguously 
able to resolve the ground state in this limit. However, the 
ambiguity is small: we show in Appendix|B]that the system at 
zero temperature is in one of only two possible phases. One 
of these is the same non-coplanar state which we find in the 
small q limit (the other is a coplanar state). This supports the 
notion that, at least up to some critical 0(1) value of q (and 
possibly for all q), the ground state evolves smoothly from the 
small q limit. 



A. Transformation to rotated frame 

From here on, we assume g <C 1. It is advantageous to 
rotate the chain spins towards the direction of the spiral mag- 
netic field 



^x,y 



l-Cx ^ X ^y 1 



+s\Ti{qx) -~cos{qx)\ 
Tlx^ \ Q - C0B{qx) - sm{qx) . (43) 

-10 / 

The rotation TZx amounts to — 7r/2 rotation about the y axis 
followed by a rotation about the z axis with angle qx. We find 
that under rotation (l43T l the magnetic field term becomes 



Hs = -2/io cos(<z/2) J2 Sx+i 



/2,y 



(44) 



Hr, transforms into 



^0 = JJ2 is:-lyS'U^.y+^O^l[SL^f2,yS'x + ^/2,y 



-Si 



x,y 
'x-l/2,y^x + l/2,y\ + '^^^ l[^x-l/2,y^x+l/2,y 

-SU/2,ySl^^/2J)■ (45) 

In the limit of small q, which we focus on, this Hamiltonian is 
conveniently split into Heisenberg one Hq with the modified 
exchange constant J — J cos q. 



Hq — J/_^^x-^.y • ^x-\ 



an effective Ising anisotropy Hi along the x axis, 

Hi^jii~cosq)j2s:.s:^^^^y 



(46) 



(47) 



^■,y 



and an effective Dzyaloshinskii-Moriya (DM) interaction H2 

H2 = Js,inq2_^ \_^x-l/2,y^x+l/2,y ~ ^l-l/2,yS x+l /2.y- 
x^y 

(48) 
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Thus, Ho = Ho + Hi+ H2 



B. Chiral SU(2) rotation 



1. Low-energy limit 

It is appropriate now to take a low-energy limit, for which 
we use the non-Abelian spin current formulation. The zeroth- 
order Heisenberg Hamiltonian in the continuum limit yields 
the fixed point term plus a backscattering correction, Hq -^ 
Ho + Hiis ■ The fixed point term, written in the Sugarawa form, 
is 

Ha = —Y^dx [JL,y ■ Jl.v + Jr,v ■ Jr,v\ , (49) 



where u = cos qu is the modified spinon velocity. To the order 
we work in this section, it is sufficient to take w w u = tt J/2. 
The backscattering correction is 



^bs = ^ dx gbsJL,y ■ JR,y^ 



(50) 



where gbs = cos qgts < 0. Again we can approximate gbs ~ 
gbs here. 

The Ising anisotropy can also be expressed in terms of cur- 
rents. One must take care since it is a composite operator. One 
obtains 



Hi = {1 — cos 5) \J / dx 



27ru 



+gbsJR^yJL,y\ ■ 



o [JL.y^L.y + JR,yJR.,y) 
(51) 



The DM (see Ref. |63| for details) and external field Hs terms 
add up to 



H2 + Hq — 



V 



dx 



d{Jly - Jly) - h{JI,^y + JIJ 



(52) 

where d = (3/7r)singJ and h = 4/io cos(g/2). We will 
consider both contributions in Eq. ( |52] ) on equal footing. For- 
mally, we consider q ^ 1 and J' /J <C 1, but with qJ/J' ar- 
bitrary. In this limit we may approximate J c:^ J, d ~ iqJ/ir 
and h ~ 2/io. Moreover, Hi can be dropped completely, since 
it represents next, q^, order anisotropy corrections to both Hq 
and the marginal backscattering term H^s- (In principle, terms 
of order q^ could be included by considering velocity shifts 
and small anisotropy corrections to the backscattering cou- 
pling gbs- However, when q is small enough, these higher- 
order corrections will not affect the outcome of the analysis in 
an essential way.) 

Finally, the relevant interchain interactions read 

H" = Y. "^-^ i^NNy ■ Ny+i + -fe£yey+i} ■ (53) 
y ■' 

The analysis in Sec.HIIlshows that coupling constants 7jv.e are 
of order {J')"^ / J^, and importantly, ^m > 0. 



An unique feature of the WZNW field theory is its emer- 
gent chiral symmetry under independent SU(2) rotations for 
the left and right moving sectors. We take advantage of this to 
remove the DM term in Eq. ([52]). Specifically, we rotate the 
right and left chiral spin currents about the y axis by opposite 
angles, +6 and —6, respectively. 



'R.y 



T^rJr. 



y 



J 



L,y 



TIlJ 



L.y, 



with 



n 



R/L 



cos{9) T 8111(6*)' 

1 
±sin(6') cos(6l) 



(54) 



(55) 



where 6 = + ata.n{d/h). Under the rotations ( |54] |. the stag- 
gered magnetization and dimerization transform to 



(56) 



Due to the chiral SU(2)flX SU(2)l symmetry , the low-energy 
Ho Hamiltonian is unaffected. Vector perturbation ( l52b sim- 
plifies to Hz, where 



ir - 


- Nr, 


Nl - 


-* coaONy + sin6'£y, 


^v - 


-* cosOsy -sinONy. 



H, 



{hf + (d)2 ^ j dx{Jly + J^,,). (57) 



Under the rotation the back-scattering term transforms into 

^''^ = E { f (JLJr.V + JL,yJty) (58) 



52 



~^~^i'^L,y'^R,y + JL.y'^R,y) + diJl^y^R,', 



53 



\'^L,yJ R,y ^ 'JL,y'^ R.y 'JL,y'^R,v '^L,y'^R,y) p 



where the couplings gi can be expressed in terms of the 
backscattering gbs 



51 



gbs 



(1 + cos 26*), 



52 



gbs 



(cos 26*- 1) 



(59) 



53 = gbs sin 26, 



.94 = gbs cos 29. 



The interchain perturbation H" is significantly affected as 
well, and now reads 

x,y 

+ sin2 e{jNeyey+i + 7e7V^iV^\i) (60) 

+ (7JV - -fe)cosesm0{eyNy^^ + A^^e^+i) 
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C. Absorption of the field h^ 

The benefit of the chiral rotation is that within Abelian 
bosonization^^i^i^ we can now absorb the magnetic field 

/i^ := J{h)^ + (J)2 by the usual shift^^ 



1 /i^ 



y^s.y 



2tt u 



■X, 



(61) 



for all y. Note that /i^ = J'^{2sqY + [iqJ/iiJ'y is 0{J') 
in the scaling limit considered here [see the discussion after 
Eq. dlOl]. 

This transforms the currents in the following way 



J. 



L,y 



J. 



j± ^±t^x 
''^L,y'^ ' 


7± ^ 7± pT'-^^ 


(62) 






(63) 



L,y 



The latter transformation explicitly embodies the linear re- 
sponse of the chain magnetization Aly = Jr y + Jl y ^'^ '^he 
field. 

The scaling dimension 1/2 fields transform according to 



Nf". 



X I iY,"; + sin 



(64) 



-x]N,^, (65) 



and Ny and N^ remain unchanged. Making this shift renders 
several terms in Eqs. ( I58I60I ) oscillatory, at scales x > u/h^. 



D. Renormallzatlon group equations 

Now we consider the effect of the various couplings. Be- 
cause of the explicit oscillatory factors introduced by the shift 
in Eq. (l6Tl i. we must consider two separate regimes of the flow. 
First, on scales shorter than the period of these oscillations, 
the oscillations themselves can be neglected, and we should 
consider all the couplings in Eqs. (15 8 1601 ). On longer scales, 
the oscillatory couplings may be dropped entirely. The reader 
may be familiar with a similar treatment of the effects of a 
field on a Heisenberg chain by Affleck and Oshikawai^ 



1. Short-scale flows 

Consider first the short-scale flows, i.e., the regime when 
—e^an < 1. This means 

71 U r^ 



no effect in this energy range, the flows remain fully SU(2) 
symmetric. They are simply 



dg. 



bs 



h 



'bs[ 



d£ 



2'KU 



(67) 



and 



d7_ 



di 


- ^^ 4..-'^^^^^' 


(68) 


d'ye 


3 




di 


- ^^+inu''^^^- 


(69) 



One can check this simple result by directly calculating the 
flow equations for gi • • • 54, and showing that the forms in 
Eqs. ( |59] l are preserved by these equations. This is the result 
of the simple argument above. 

We note that it is sufficient to work only to linear order in 
the relevant couplings jn, 7^, since their initial values are of 
order (J')*, and therefore remain small over this range of i 
(they increase only by a factor of e^ ~ J/ J')- To this order, 
they do not feed back into the flow of g^s- The usual solution 
to Eq. ( |67] i obtains: 



9bsii) 



9bs 



a-§^ey 



(70) 



Since gi,s < 0, it becomes small under renormalization, and 
specifically of order u/£ for i ^ 1. Inserting this into the 
remaining equations and solving gives 



lN{i) 



p\ 1/2 
, 9bs*- \ i /„N 

1 - - — I e 7Ar(0), 



2t:u I 



(71a) 



-3/2 



le{t) - (l-f^) e^7e(0). (71b) 



Evaluating this at 
conditions gives 



= (.* and using Eqs. JTSI l for the initial 



7e 



1 



2'KU 

\gbs\HJ/J' 



.P 



^(\g,s\HJ/J')\^'\j'f 
27r4 I, 

3 

^4^ 



2'KU 



(J' 



J2 



(72a) 
(72b) 



Note that these couplings indeed remain small at this scale. 
Furthermore, the staggered magnetization coupUng 7jv is al- 
ready parametrically enhanced over the dimerization coupling 
7^ , by a factor of In" ( J/ J' ) . 



< ^ < /* = In ■ 



h^ao 



MJ/f) 



(66) 



We neglect completely the effect of the oscillatory factors in- 
duced in Hts and H". Although the form in Eq. ( fSSl l is com- 
plicated, the flows remain simple. This is because Eqs. (|58l59t 
are obtained from the chiral SU(2) rotation which is a symme- 
try of the fixed point Hamiltonian. Thus, since the field /i^ has 



2. Long-scale flows 

Now we consider the renormalization on scales longer than 
the period of the oscillations induced by the field shift. Here 
the SU(2) symmetry is truly broken, and the RG deviates from 
the simple one above. We drop all oscillating terms (this in- 
cludes 51, (73 and several of the terms in H"), so that the re- 
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maining perturbations to the bare Hamiltonian Hq are 



a:,y 



^'"^ " Hiy^-^i.y^i.y 



■h,yJR,y) 



94:JL,yJR,y 



(73) 



x,V 



+7+(7V^^7Vy\l+£,yey + l)}, 

where we have introduced the new coupUng constants 
92,94, IN^ , jNy , 7+- They should be matched at £ = £* to the 
couplings from the short-scale flows, defined in Eqs. ( I60I58I I. 
which implies 



hin 


= -gbs(r)sin2 0, 


un 


= ,9bs(r)cos20, 


7iv-(r) 


= lN{t), 


lm{t) 


= cos2 07jv(r)+sin2 07ar), (74) 


1+in 


= cos2 07,(r) + (1 + sin^ 0)7jv(r)] /2. 



We now compute the RG equations for the perturbations (l73l l 
to the bare Hamiltonian iJo using well-established OPEs for 
the non-Abelian spin currents (see for example Re f .140). After 
some lengthy calculations we find 



dg2 
di 

d£ 

d^Nv 

di 

dn+_ 
di 



9492 



dg4 
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= 1 
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At:u 

1 
At:u 

1 



94 



94 



1 

2t:u 

1 
2t:u" 



92 I IN-' 



92 iNy , 



(75) 



94 7+- 



It is important to understand how these equations lead to an 
instability. The equations for g2 , g4 are decoupled, and can 
be solved separately. They have the standard form found in 
the Kosterlitz-Thouless (KT) analysis42i One recalls that the 
quantity 



:^2 



Y = g'2-~9i 



(76) 



is a constant of the motion. The flows are unstable provided 

g4{i*) < \g2{i*)\, which is always satisfied except when 6 = 
7r/2 exactly, at which point it becomes an equality. That is 
for 9 G [0, 7r/2[, the trajectories tend to 94 -^ —00 and 32 ^^ 
+00. We fix Y by the initial conditions. 



Y^\gUn\' 



(sin"* 9 - cos^ 26 



(77) 



Hence, Y is negative for 9 G [0, acos a/2/3] and positive for 
9 g [acos \/2/3, 7r/2]. Writing .92 = ^ + .94 we can solve 
the KT equations for 94. When Y > 0, i.e., in the crossover 
regime of the KT flow, we have 



atan 



94{ 



atan 



gbs{i*)cos29 



Y 



2tiu^ 



r). 

(78) 



The coupling g^ clearly diverges when the right-hand side of 
this equation reaches 7r/2 plus an integer times tt. The "time" 
id of this divergence is 



27ru 7r/2-atan(i^) 



(79) 



where we define T = \/Y /\gb,{i*)\. Using i* - ln(J/J') 
and gbs{i*) ^ 2TTu/i* we obtain 



id = ln(J/J') 



7r/2 — atan | 



i2e\ 



(80) 



One can check that id/ h\{J / J') increases monotonically 
from 4 when 9 = acos^2/3 to infinity as, 9 ^ -k 12. Sim- 
ilarly, when y < (strong-coupling regime) 34 (^) diverges at 
the length scale 



id = HJ/j') 



1 + atanh 



T 



cos 26* 



(1/T) 



(81) 



where T = ^/^ /\gi,s{i*)\. In this regime ^rf/ In (J/ J') 
takes the value 4 at = acos ■\/2/3, increases monotonically 
with decreasing 9, and diverges at ^ = 0. 

With the values ati = i* given by Eqs. (|72]|74]i the relevant 
couplings 7jYx, 77VK> ^nd 7+ become of order 1 at the length 
scale io ~ 41n(J/J') [so that ^ - T - 31n(J/J')], which 
is always smaller than the scale id- The most relevant of these 
turns out to be 7Arx, as can be seen, e.g., by examining the 
ratio 



7iV-ffl ^ iN-ji*) 
lNy{i) lNy{i*) 



exp 



1 



■nu Ji 



dxg2{x) 



(82) 



where (remember that g^s < 0) 



hin = \gUi*)\sin^O 



\9bS*)\d' 
d? + K^ 



(83) 



From this clearly 7Ari becomes parametrically larger than ^^y 
under renormalization (for q ^ 0). Similar analysis shows 
that 7Arx is also enhanced relative to 7+. Thus, for 9 > 0, 
we find that the staggered magnetization 7Arx = 77V- domi- 
nates. Hence, the spins align antiferromagnetically along the 
X direction, S^, in the rotated ("comoving") coordinate frame. 



3. Ordering pattern of chain spins 

Let us now infer what this means in terms of the orig- 
inal spins, in the fixed coordinate frame. It is necessary 
to trace back the transformations of the spin operators in 
Eqs. ( 143 15416 il l. This is straightforward but tedious. We will 
not give a general expression of the relation of the microscopic 
spins to the continuum operators after the final transformation, 
which is not illuminating. Instead, we give the result for the 
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expectation value of the spin operators given that, as argued 
above, in the rotated variables the ordering is very simple: 



{N^) - 


= M{-iy6'''\ 


(84) 


(Jb) 


= (J£) = (e) = 0. 


(85) 



Here A/ ^ represents the spontaneous moment, and will be- 
come the staggered magnetization. The (—1)^ factor obtains 
because 7Arx > 0. 

Now we relate the spin operators as described above to the 
continuum fields: 



(St 



x+-!j,yi 



^IhJ 



2ttu 



cos 9 cos gx, 



cos a sin qx, 



{s: 



+^.y' 



2ttu 



(86) 

(87) 
(88) 



Thus, indeed the dominating 7jv^ -term orders the chain spins 
antiferromagnetically along the z axis perpendicular to the 
x-y s-spiral plane. The non-zero components of the spins 
within the x-y plane are induced by the local field arising 
from the ordered interstitial moments. Note that they are of 
0{J'/J) which is much larger than M -- {J'/J^- Thus 
the static moments on the chains are predominantly in the the 
plane of the spiral, with a smaller staggered component in the 
perpendicular (z) direction (see Fig.|2l). 

The case of ferromagnetic order amongst the interchain 
spins is obtained by setting q = 0. This implies d = so 
that 32 = 0, see ( [83] l. The chains are subjected to the uniform 
magnetic field h only. This leads to 9 = and, as a result, 
symmetry in the N^ — N^ plane: '^n<^ /iNy = 1- The chain 
spins S order nearly collinearly with a ferromagnetic compo- 
nent of order 0{J' / J) and smaller antiferromagnetic compo- 
nent of order 0[{J' / JY] perpendicular to the ferromagnetic 
moments. Note that in this case, because the interstitial spins 
and the predominant moment of the chain spins are ferromag- 
netically ordered and hence collinear, the full magnetic order 
is actually coplanar. For instance, if the ferromagnetic mo- 
ments are aligned in the x direction, the antiferromagnetic 
component of the chain spins will be aligned along some axis 
e in the y-z plane, and all the spins are contained in the plane 
spanned by x and e. Furthermore, this state is ferrimagnetic, 
i.e., has a macroscopic net moment, since the moments of the 
ferromagnetically ordered interstitial spins are unequal to the 
opposing in-plane component of the chain spin moments. 



VII. SUMMARY AND DISCUSSION 

In this work we have analyzed the ground state phase of 
the quantum Heisenberg antiferromagnet on the kagome lat- 
tice with spatially anisotropic exchange. We have studied this 
problem in the quasi- ID limit, where the lattice is broken up 
into antiferromagnetic spin- 1/2 chains that are weakly inter- 
acting via intermediate "dangling" spins s (see Fig. ^. This 
limit lends itself to a perturbative RG analysis in the weak ex- 
change interaction J' using bosonization and current algebra 



techniques. We find that there is a natural separation of en- 
ergy scales: the intermediate spins order at an energy scale of 
order {J'Y / J, at which the chains are not influenced by the 
interactions among themselves and with the interstitial spins. 
The low-energy behavior of the chains, on the other hand, 
is only modified at 0[(J')^], as geometric frustration pre- 
vents the generation of relevant interchain interactions at the 
larger scale {J'Y / J. We have used perturbative and numeri- 
cal methods to determine the effective interactions Ji among 
the interstitial spins s, which arise at 0[(J')^]. It turns out 
that the spins s order in a coplanar cycloidal spiral with wave 
vector q parallel to the chain direction. The magnitude of the 
wave vector is presumably rather small (if not vanishing). It 
depends sensitively on the strength of further neighbor inter- 
actions among the interstitial spins, which cannot be reliably 
determined from our approach. The ordered interstitial mo- 
ments induce a spiral order of the chain spins of 0{J' / J). 
Besides this, the chain spins exhibit a small antiferromagnetic 
component of 0[(,/'/,/)^] that points along the axis perpen- 
dicular to the spiral plane. This non-coplanar ground state of 
the spatially anisotropic kagome antiferromagnet is illustrated 
inFig.|2] 

It is interesting to compare to recent results for this lattice in 
the spatially anisotropic limit. Wang et alM found a coplanar 
ferrimagnetic chiraUty stripe order using a semiclassical anal- 
ysis (see also Ref. l37h . In this state, the interstitial spins are 
ferromagnetically ordered, and the chain spins are ordered in 
an antiferromagnetic fashion, nearly collinearly along an axis 
perpendicular to the interstitial moments, but canted slightly 
in that direction. The ordering of the interstitial spins is very 
close to our findings; i.e., it corresponds to the special case 
q = 0, which as we have described cannot be excluded by our 
calculations. However, even in this case the ordering pattern 
of the chain spins is quite different from that in Ref. [381 insofar 
as we find that the chain spins have a predominant /erroma^- 
netic component antiparallel to the interstitial spins (and only 
a considerably smaller antiferromagnetic component perpen- 
dicular to the interstitial spins), while Wang et ali^ obtained 
a predominantly antiferromagnetic ordering among the chain 
spins. 

Yavors'kii et al. in Ref. [33, on the other hand, used a large- 
N expansion applied to the Sp(-/V)-symmetric generalization 
of the model. In the Umit J' <^ J they found that the chains 
are completely decoupled, and the interstitial spins show some 
(short range) spin-spin correlation that is compatible with a 
spiral ordering pattern. While the mean field treatment of this 
large-A^ approach seems to miss the predominant spiral or- 
dering of the chain spins, the spiral ordering of the interstitial 
spins is in agreement with the findings of this work. 

Numerical studies of the spatially anisotropic kagome 
model should be very helpful in establishing the range of 
spatial anisotropy of exchange interaction where the non- 
coplanar ordered state found in this work represents the 
ground state of the system. We hope that our work will in- 
spire further investigations of this interesting problem. 
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say. These local contributions read 

-T J dTidr2s;;+i(Ti)sj;_^_i(T2) 



(Al) 



+ 27172 [M-{t,) + Af,%i(Tl)] d, [7V^^T2) + N'y+,iT2)] 



APPENDIX A: RG CALCULATIONS 



In this Appendix we present a detailed derivation of flow 
equations (fT4l) and ( fTsl ) using standard OPE relations for the 
continuum fields. 



1. Derivation of Eqs. fl4t 

The first order terms in the RG equations ( fT4] i originate 
from the rescaling of the time and space coordinates and the 
redefinition of the fields [see Eqs. (fTZt l. As already discussed 
in the main text, the interaction terms Vi and V2, Eqs. (fT3T l. are 
local in space, which means that their scaling dimensions have 
to be compared with 1, the dimensionality of time r. Intersti- 
tial spins have scaling dimension 0, while M (N) fields are of 
scaling dimension 1 (1/2). Consequently, the scaling dimen- 
sion of Vi is 1, while that of V2 is 3/2, from which follows that 
the RG equation for the coupling 71 does not contain a linear 
term, while that for 72 starts with (1 — 8/2)72 = —72/2, in 
agreement with the first line of Eqs. (fl4] l. 

The second and third lines of Eqs. ( fT4] l describe the flow 
of the fluctuation-generated interchain couplings 14h, Eq. (O, 

and Vpjj , Eq. (|9]l, which operate in two-dimensional space- 
time. As a result, the scaling dimensions of the effective in- 
terchain interactions have to be compared with 2. The first 
order terms in the second line of Eqs. ( fT4l i are then a direct 
consequence of the fact that the scaling dimension of the 7on- 
interaction is 3, while that of the 7j\f -interaction is 2. Sim- 
ilarly, positive linear terms in the last line of Eqs. ( fT4l i fol- 
low from the strong relevance of the interchain couplings V"ch, 
Eq. dS), which have scaling dimension 1. 

The second-order corrections to flow equations (fT4] i are de- 
rived from contracting terms in perturbation expansion ( fTTT i. 
In order to obtain the one-loop corrections in Eqs. ( fT4l i. we 
need to consider contributions to the second order term in 
Eq. ( fTTl i that either yield a renormalization of the couplings 
ji or generate new interchain interactions. We begin by se- 
lecting from these contributions terms that contain a product 



where we have suppressed the x coordinate for brevity. Even 
though the interchain spins s _|_j. have no dynamics of their 

own at this level, they must be time-ordered as follows^ 



Ts" 



= 0r,-r,S^iTl)s\T2) + 9 r, - r^ -s" {T2) S^ (t,) 

5°-'' i 



H 



^abc ^c 



^^W, 



(A2) 



where 9t is the step function, t = (ti + T2)/2 is the center- 
of-mass time, and i := ti — T2 is the relative time. 

The off-diagonal term in Eq. ( IA2b (which is proportional 
to e'^''^) is responsible for the renormalization of the Kondo- 
like couplings 71 and 72. As an example we consider the 
renormalization of 71, which comes from the second line in 
Eq. (lAU . Separating slow and fast degrees of freedom we 
can apply OPE (|4|i (and a similar expression for the left cur- 
rents) to the product of two spin currents at nearby points [i.e.. 



(IA2Ti leads to 



7? 



ie"'"^A/^(T)/(27ri)]. Combining this with 



4ttu 



dTslMT)M':{T) 



v+ 



a<\t\<ba 



dt 



(A3) 



The integral with |i| > ba does not contribute to the renor- 
malization. The one-loop correction to the flow equation for 
7i can now be read off (IA3l l as ex dl^l/iru, which gives us 
the first equation in (fl4] i. 

The renormalization of 72 is computed in a similar way, 
one only needs to realize that it comes from the third line in 
Eq. dATJ. Fusing M^ with d,^N^ via the OPE (see Ref.H^for 
more details) 

M'^{Ti)d,Nl{T2) = lim M!^{x' ,Ti)d^N'^{x,T2) 



-(5"%(x,r) le'^'^^d^N'yix.T) 



27r{ut + aoCFty 2TT{ut + aoat) 



(A4) 



of two intermediate spins from the same site, {2x - 



,y- 



leads to the following one-loop correction of the flow equation 
for 72: 5j2 oc (i£7i72/7r. It is useful to note that rescaling of 
space and time does not affect quadratic terms, as each of them 
is explicitly proportional to the RG step di, which comes from 
the shell integration of relative coordinates. 

The second-order corrections to the flow equations for the 
interchain couplings jm and -fgN follow from the diagonal 
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term in Eq. ( IA2I ) (which is proportional to (5°^). For example, 
applying relation ( IA2I ) to the second line of Eq. ( lAlb produces 



^JdrM^{r)M^^,{. 



dt. 



Q<|t|<6Q 



To generate from this the ^m term in Eq. (|9]l one needs to 
sum all local contributions like the one above using 'Ylx ■■■ ^ 
/ dx/{2ao)... , as appropriate for the kagome geometry. As 
a result one finds that 7^/ ^ jfdl. Similarly, the other inter- 
chain coupling, jgN in Eq. (|9]l, can be derived starting from 
the last line in Eq. (lAlb . 

Finally, we turn to the relevant interchain interactions 7„ 
and 7^ which are generated by fusing the marginal interaction 
7j\f and the irrelevant interaction 7ojv Details of this pro- 
cedure have been described previously in Refs. I39p0i Here 
we would only like to mention that the RG scheme that we 
have adopted here (i.e., integrating the one-loop x integrals 
over the entire space of relative x coordinates while restrict- 



ing the relative time integral to the shell, a < It 
which is different from the RG scheme of Refs. |39 



< ba), 
I40l does 



not modify the outcome of the calculation in Refs. I39p0l in 
any significant way. Namely, upon fusing A/"(a;i,ri) with 
dxNy{x2,T2) on chain y (and, similarly, on chain y ± 1), 
one arrives at the following integral over the relative coor- 
dinate x = Xi — a;2 and over relative time t = ti — T2: 
I - lZd^L<w<t.dtix' + er' -{h- l)a-\ This 
explains the structure of the quadratic terms in the last line of 
Eqs. (fl4] i. 

Additionally, we note that the third line in Eq. jAlb in com- 
bination with the first term in Eq. ( IA4b results in a strongly 
relevant contribution (scaling dimension 1/2) 



-67172 
■naou 



dTey{x, t). 



(A5) 



In the two-dimensional version of the spatially anisotropic 
kagome lattice this term cancels out, since the summation 
over local bow-tie crossings (that is, over x) brings in the 
factor (—1)^ [which originates from the staggering factor in 
Eq. (I2ali1. resulting in /dTdx(— l)^£j^(a;, r) -^ 0. This is 
how symmetry (l3cT l manifests itself. 

In contrast, for the kagome strip of extension one in the 
y direction, the staggering factor (—1)'^ does not appear, 
since the bow-tie crossings are separated by two lattice spac- 
ings. As a consequence, the expression (IA5b turns into 
/ dTdxey{x, r), which implies spontaneous dimerization of 
the kagome strip, an ordering pattern that does not reduce 
the translational symmetry. This strongly relevant term was 
missed in previous analytical work^~. Numerical studies, on 
the other hand, did find dimerized ground states^. 



2. Derivation of Eq. QD 

To derive flow equation ( fTSl l we need to select from the sec- 
ond order term in Eq. ( fTTT i contributions that contain products 
of two different intermediate spins with the same y coordinate, 
2/ + ^, say. Among these, the most important contributions are 



those, that involve products of uniform or staggered magneti- 
zations from the same chain 

+2x^T^ /rfndr2 4^^^.(Ti)4^+i(r2) 

x[^lAr{2xi,Ti)M\2x2,T2) (A6) 

+-fldxN'' [2x1, Ti)dxN'' {2x2, T2) 

where we have suppressed the y coordinate for brevity. The 
factor of 2 in the first line arises because there is an equal 
contribution from both the y and the {y + I) chains. Since 
the spins Sj^-^i at different sites commute and time-ordering 
of the continuum fields M and N is automatic, we can dis- 
regard the operator T in the above expression, provided we 
exclude the case xi = X2, which was treated in the previous 
subsection. 

By splitting the integrals of Eq. ( IA6I ) into slow and fast de- 
grees, we can use the OPEs, Eqs. (|4]i, to fuse the product of 
two continuum fields at nearby points. In this way, we derive 
from Eq. (IA6b the one-loop renormalization to the first term 
of the interaction i^A, Eq. ( fTOb 

E / ^^4., + i4..+i (if III + iIIm + iPn) ,(A7a) 



with the integrals 
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dt- 



l/(87r2) 



a<\t\<ba [ut ± i{xi — X2) + a^at 



Q<|t|<6Q 



dtdx^dx^ 



^JuH^ + A{xi -X2) 



2, (A7b) 
<A7c) 



where Cat w (27r)~^/2 is the amplitude of the {N°-N°-) cor- 
relator. The integral over the infinitesimal interval [a, ha] in 
Eqs. ( |A7bb and (|A7cb amounts to replacing t with a^. The 
subsequent rescaling turns the cut-off ai back into the mi- 
croscopic cut-off a, while x —> x^ = xe~^ . Simplifying 
/j^j and explicitly taking derivatives in /^v leads 



t-M 



'A/ 



to Eqs. ( fTSl l. It is worth noting that the two contributions Im 
and /tv are of opposite signs, resulting in a fast decay of inter- 
actions between interchain spins s. 



APPENDIX B: CHAIN ORDERING IN THE LIMIT OF A 
LARGE SPIRAL WAVE VECTOR (q ~ 0[1]) 

In this Appendix, we briefly discuss how the chain spins 
might order under the perturbation of a spiral magnetic field 
with a large wave vector q = ((7,0), where q ^ 0\\\. In 
this case the transformation to a rotated frame as done in 
Section rVlAI is not useful, as the generated interaction terms 
come with couplings of the order of the bare exchange J, see 
Eqs. Wh and ( |48] |. Instead, we should remain in the original, 
non-rotated basis, which leaves the dominant 0( J) interac- 
tions in their simplest form. Moreover, the rapid oscillation of 
the spiral field, which is highly incommensurate with the "nat- 
ural" wave vectors and -n of correlations of the Heisenberg 
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chains, ensures that its effects rapidly average out to leading 
order However, to 0{h^) we may expect non-oscillating in- 
teractions to be generated. In principle these can be obtained 
by a perturbative analysis expanding in powers of Hs ( l40l i. 

As usual, symmetry analysis is helpful in figuring out the 
type of terms that can be expected from such a calculation. 
There are two important symmetries. First, the spiral field 
term breaks both spin-rotational and translational symmetries, 
but is invariant under a translation x -^ x + \ followed 
by a simultaneous spin rotation by the angle q about the z 
axis. Spiral state ( [39] l is also invariant under spatial inversion 
{x -^ —x) followed by a change of sign for the y and z com- 
ponents of the spins (5''''^ -^ — S*^'^). Taking into account 
these constraints, the only possible marginal or relevant terms 
which may be generated in a single chain are 



xy 



Jdx {d^ir^j, - J^l) + Sg,J^j,J^L} ■ (Bl) 



We neglect here terms that are already present without the spi- 
ral field, and those that couple different chains as these cou- 
pling constants are necessarily smaller by at least one power 
of J', since the latter cannot be generated without some bare 
inter-chain interactions. 

A naive current algebra calculation using the continuum ap- 
proximation, Eq. ( |2al i for the spin operators, indeed produces 
precisely these terms at 0{h^). However, the resulting ex- 
plicit expressions are not reliable, as they are dominated by 
short distances of order q^"^, in which the continuum limit is 
inappropriate. A lattice scale calculation, similar in spirit to 
that in Sec.|IV](yet more closely to that in Ref42), is required. 
Unfortunately, it turns out that to do so requires detailed infor- 
mation on the lattice scale properties of certain four-spin cor- 
relation functions of a Heisenberg chain. These data are not 
available to our knowledge. Therefore, we must rely upon the 
symmetry considerations alone, assuming no particular signs 



or magnitudes for d^ and Sg^ above, apart from the fact that 
both are expected to be of 0{h? / J). 

Fortunately, this does not result in significant ambiguity. 
This is largely because the nominally "relevant" DM correc- 
tion dz has trivial effects. Similar to the magnetic field in 
( |6Tl i it can be removed by a simple shift. Unlike the case in 
Sec. 1111 B1 however, the shift does not affect the backscattering 
H\,s, which is written in terms of the field ip, which is dual to 6. 
The only effect of the shift is to change the ordering momen- 
tum, if any, of the N^'^ components. It does not determine 
the nature of the ordering instability. 

The anisotropy term 5gz is important, as it tips the balance 
of competition between different components of TV fields in 
the interchain Hamiltonian H", Eq. ( |53] |. With the help of an 
OPE-based calculation similar to the one that led to Eq. ( l76b 
we find 



d_ "IN' _ Sg^ 
di 77v^ At:u 



(B2) 



This tells us that the type of AT-order is determined by the 
sign of the generated Sgz- Positive 6gz favors N^ compo- 
nents, leading to the non-coplanar ordering pattern found in 
Section [VID 31 in the small-g limit. Negative 6gz, on the 
other hand, would prefer N^-^ components, without breaking 
the symmetry between them. Such a state is clearly co-planar, 
and different from the one found in Section 1111 B I Note that 
the DM term dz will affect the ordering wave vector of this 
state but not the non-coplanar one. Which of these two sit- 
uations is obtained cannot be discriminated by our analysis, 
since the sign of Sgz is not determined. However, it is rather 
natural to expect that as q is reduced from 0(1) values, one 
should observe behavior consistent with the small q analysis. 
This would suggest that Sgz > for a non-vanishing range 
of q greater than zero, and indeed it is possible that this is the 
case for all q. 
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